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ABSTRACT 

Dwarf spheroidal galaxies (dSphs) are extremely gas poor, dark matter-dominated galaxies, 
which make them ideal laboratories to test the predictions of the Cold Dark Matter (CDM) 
model. We argue that the removal of a small baryonic component from the central regions 
of forming dSphs may substantially reduce their central dark matter density. Thus it may 
play an important role in alleviating one of the problems of the CDM model related with 
the structure of relatively massive satellite galaxies of the Milky Way. Traditionally, colli- 
sionless cosmological iV-body simulations are used when confronting theoretical predictions 
with observations. However, these simulations assume that the baryon fraction everywhere in 
the Universe is equal to the cosmic mean, an assumption which can be incorrect for dSphs. 
We point out that the combination of (i) the lower baryon fraction in dSphs compared to the 
cosmic mean and (ii) the concentration of baryons in the inner part of the Milky Way halo 
can go a long way towards explaining the observed structure (circular velocity profiles) of 
dSphs. We perform controlled numerical simulations that mimic the effects of baryons. From 
these we find that the blowing away of baryons by ram pressure, when the dwarfs fall into 
larger galaxies, lowers the circular velocity profile of the satellite. In the likely scenario that 
baryons are blown away after the central ^200-500 pc region of the galaxies was formed, 
we show that the dark matter density in the central regions of dSphs declines by a factor of 
(1 — /b) 4 ~ 0.5, where /b is the cosmological fraction of baryons. In addition, the enhanced 
baryonic mass in the central regions of the parent galaxy generates tidal forces, which are 
larger than those experienced by subhaloes in traditional TV -body simulations. Increased tidal 
forces substantially alter circular velocity profiles for satellites that come as close as 50 kpc. 
We show that these two effects are strong enough to bring the observed structure of dSphs 
of the Milky Way into agreement with the predictions of the subhaloes in CDM simulations, 
regardless of the details of the baryonic processes. 

Keywords: cosmology: theory -dark matter -galaxies: haloes - galaxies: dwarfs - methods: 
N-body simulations. 



1 INTRODUCTION 

According to the concordance cosmological model, a large frac- 
tion of the matter in the Universe ought to be collisionless and 
dark. This Cold Dark Matter (CDM) model is supported by the ob- 
served large-scale structure in the Universe, as traced by galaxies, 
and the emergence of this structure from the level of initial density 
fluctuations observed in the cosmic microwave background (e.g., 
iKomatsu et alj l201~ll ). However, some inconsistencies between a 
number of observations and the predictions of the CDM model on 
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small (galactic) scales have understandably produced concern and 
attracted a great deal of attention from both observers and theorists. 

The first of the small-scale problems for the CDM model 
is the so-called "m issing satellites problem" jKlvpin et al.ll 1 999bl : 
iMoore et al.lll999l) . Numerical simulations of the CDM theory pre- 
dict a large number of low mass substructures to be present in 
a dark matter halo comparable in size to that of the Milky Way 
(MW). However, the number of subhaloes is much larger than the 
number of satellite galaxies found around MW and in the Local 
Group, even after accounting for the extreme low surface brightness 
dwarf galaxies th at have been recently found using the Sloan Dig- 
ital Sky Survey jTollerud et al.ll2008fc IKqposov et ai]|2008l. 120091 : 
lBullocketalJl2010h . 
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These observations of the abundance of satellite galax- 
ies can be reconciled within the context of CDM by in- 
voking galaxy formation physics. Star formation in low 
mass subhaloes can be supp ressed by the heating of cold 



gas due to photoionization ( Quinn, Katz, & Efstathiou 
Navarro & Steinmetzll 1997k iBullock. Kravtsov. & Weinberg! 



1996 



2001 



Somervillej 120021: iBenson et alj |2002| ; lOkamoto. Gao. & Theunsl 
20081 : lOkamoto & Frenkl l2009h and in large r subhaloes through 
energetic feedback from sup e rnovae (e . g.. iDekel & SilkJ [l98i 



Navarro. Eke. & Frenkl Il996at iKravtsovl BOlOl IPefiarrubia et al j 



2012b . Thus the solution to the missing satellites problem may be 



in the complex physics of galaxy formation. However, it is worth 
noting that models with realistic and physically-motivated prescrip- 
tions remain challenging to accurately simulate due to the high 
numerical resolution required to capture the relevant physical pro- 
cesses. 

A related issue is the abundance of relatively massive 
satellites, the so-called "too big to fail" problem - does the 
CDM theory predict too many satellites with large circular 
velocities (n c i rc )? It has been argued that the theory pre- 
dicts the correct abund ance of very large satellites, such as 
the Magellanic Clouds jBusha etalj l2011a!lbT). and dwarf field 
galaxies with ii c j rc > 



50 kms -1 jTikhonov & Klypirj 120091 ; 



iTruiillo-Gomez et alj 1201 ll) . However, the situation is different 



for satellites with 30 km s 



< v c 



< 50 km s . Compar- 



ing the circular velocity profiles of subhaloes in the Aquarius 
TV-body simulations with the circular velocities measured at the 
half-light ra dii of the bright Local Group dwarf spheroidal galax- 
ies (dSphs), lBovlan-Kolchin. Bullock. & Kaplinghaj t201ll . l2012h 
find that there are too many subhaloes. For example, haloes with 
Mvir ~ 1.4 x 10 12 Mq have ~25 satellites with infall circular 
vel ocities ti c j rc > 30 kms" 1 within 300 kpc of the halo's cen- 
tre jBovlan-Kolchin et alj|2o"l2h . In a given halo, at least seven of 
these satellites do not have observed MW counterparts, so the sim- 
ulations contain too many very compact, large subhaloes. 

There are a couple of solutions, within the context of the CDM 
model, that can explain some of these inconsistencies. First and 
foremost, the CDM cosmology adopted for the Aquarius haloes 
uses as = 0.9 instead of the lower concordance value of as = 0.8. 
Subhaloes in lower erg models have lower abundances and their 
density profiles are less concentrated, such th at they are more sus- 
ceptible to tidal effects jDi Cintio et al.|[2oTTh . Secondly, the abun- 
dance of substructures above a certain circular velocity scales lin- 
early with the host halo mass. Lowering the mass of the MW by ~ 
20 per cent (i.e., from M vil ~ 10 12 M to M vil sC 8 x 10 11 Mq) 
can reduce the abundance of the massive, dense subhaloes to a 
number low enough so that the absence of such satel lites around 
the M W can simply be attributed to a statistical fluke jWang et"al] 
120121) . 

The other possibility is that baryonic processes, such as stellar 
feedback and photoionization, in combination with the subsequent 
tidal stripping of their outer regions, impacts the inner structure 
of massive subhaloes. This could result in the generation of cores 
in their central densit y profiles and simultaneously solve the mas- 
sive failures problem dMashchenk o. Couchman, & Wadslev 2006 ; 



Mashchenko. Wadsley. & Couchmanl 2008 ; Governato et al 201C 
ltt 



Penarrubia et al. 12 010; Di Cintio et allfe oi 1: Pontzen & Governatol 
2012l;lGovernato et aLkoLllPari-v et al.ll2012l;lzolotov et al.l2012l; 
Brooks & Zolotovl 120121 : iBrooks et alj |2012| ; IPefiarrubia et alj 
2012b . However, capturing the baryonic processes accurately re- 
quires numerical simulations with hydrodynamics to be run at res- 
olutions of tens of parsecs, which are difficult to achieve. Ques- 



tions also remain about the sub-grid models implemented to mimic 
baryonic processes and the magnitude of the resulting feedback 
dCeverino & Klypin 2009; Cohii et alJl201Ch , 

If the discrepancies between predictions of the CDM model 
and observations listed above are not resolved with realistic bary- 
onic effects, then these observations will provide a unique op- 
portunity to probe the nature of dark matter. Alternative sce- 
narios, such as warm dark matter or self-interacting dark mat- 
ter, have been proposed to explain the diffe rences between the 
CDM model predictions and obs ervations jLovell et~aT 
IVogelsberger. Zavala. & Loebll2012h . 

In this paper, we investigate modifications to the CDM ex- 
pectations induced by the presence of baryons. We focus on well- 
established and undisputed facts about the behaviour of baryons in 
the MW, namely (i) that the fraction of mass in baryons within sub- 
haloes that host dSphs is substantially below the cosmic mean and 
(ii) that there is a baryonic concentration within the central parts of 
the MW halo. 

Cosmological TV-body simulations of large-scale structure are 
initialised taking into account the effect of baryons on matter 
fluctuations. However, while evolving this simulation forward in 
time, the simplifying assumption that all matter, including bary- 
onic matter, is collisionless, is commonly used. This implies that 
the baryonic distribution follows the matter distribution in an un- 
biased manner. Although this assumption works well if one wants 
to calculate large-scale properties of the mass distribution (scales 
larger than typical haloes), the predictions of such simulations on 
smaller scales must be inteipreted with caution. Gas poor objects, 
like dSphs, should have larger dark matter masses than baryonic 
masses, while the central parts of the MW should have more bary- 
onic mass than it is assigned in simulations. 

This simplifying assumption leads to a twofold effect, caus- 
ing the dSphs to be more resistant to tidal stripping in conventional 
TV-body simulations than in reality. The density of dSphs in the in- 
ner regions is larger in traditional TV-body simulations than what it 
should be, given the loss of baryonic material. At the same time, 
the density in the central regions of the parent halo is underesti- 
mated. In reality, baryons lose energy through radiation and fall to 
the centre of the parent halo which makes their interiors more mas- 
sive. In TV-body simulations, this is neglected because the baryons 
are assumed to be collisionless, similar to dark matter. This implies 
that the tidal forces that act on the subhaloes that host the dSphs 
are weaker in these simulations than in reality. This causes massive 
subhaloes to preserve their structure and high circular velocities, 
allowing them to survive in simulations down to redshift zero. 

The paper is organized as follows. We discuss different as- 
pects of dSph galaxies and possible scenarios for their formation in 
Section [2] In Section [3] we present the details of our simulations 
and analytical methods. The effects of baryon removal for isolated 
galaxies are studied in Section|4] In Section|5] we present results of 
tidal stripping models, and conclusions are given in Section[6] 



2 FORMATION AND EVOLUTION OF DWARF 
SPHEROIDALS 

The origin of dSphs is, as yet, unknown with current formation 
scenarios constrained by simulations and by observations of dwarfs 
in the Local Group. In order to address why dSphs have a lower 
baryonic fraction than the cosmic mean, we must examine their 
formation and evolution. This information will guide us when we 
set the parameters of dwarf galaxies in our simulations. 
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The radial distribution of satellites may provide valuable in- 
formation on the evolutionary history of dSphs. It is known that 
95 per cent of dSphs are located within distances of 250 kpc 
from the centres of known groups (e.g.. [ Karachentsevl 120051 ; 
iKarachentsev. Karachentseva. & Sharinal I^OolT In other words, 
there are hardly any isolated dSphs. Even inside the Local Group 
the dSphs are preferentially found in the inner regions around the 
MW and M31. These results imply that there are unique environ- 
mental forces altering the majority of dSph morphologies. 

A possible mechanism for forming dwarf galaxies is through 
the smallest "baryonic building blocks," which are 3u fluctua- 
tions at z ~ 10, on the order of 10 s Mq. Haloes with M < 
2 x 10 7 M are unable to retain their gas when SNe feed- 
back and reionization are introduced, but the remai ning structures 
with M > 2 x 10 7 Mq do retain their gas (see [Peebles! 1 1984 



Ricotti & GnedirJ 120051; jGnedin & Kravtsovl l200d: iMoore et al] 
20061 ; iRead. Pontzen. & Vielll2006a! ; IPefiarrubia et alj|20ld, fora 



discussion of the critical mass scale and star formation thresholds). 
When these smallest luminous structures remain far from poten- 
tial perturbers, they retain their gas and form stars, matching the 
field dwarf irregular (dlrr) population. The satellites with pericen- 
tric passages close to their host are subject to ram pressure strip- 
ping, which removes the gas and inhibits star formation. 

A formation scenario that naturally rec overs the radial dis - 
tribution of dSphs is termed "tidal stirring" dMaver et al.ll200l"h . 
The tidal stirring mechanism relies on the existence of field dlrr 
galaxies that begin discy and gas rich and enter the gravitational 
potential of a massive host before z ~ 1. Most of the dwarf's 
gas content is lost via ram pressure stripping, hindering future 
star f ormation (e.g.,[Eskridge 1988; Grebel. Gallagher. & Harbeckl 
120031) . The morphology of the dwarf is changed by tidal forces; 
the disc transforms into a spheroidal component due to bar in- 
stabil i ties created by tidal shocks (see iMaver et al. 1200 it iMaveJ 
20051; iMaver et al.l |2007|; iKlimentowski et alj 12001 120091 . 1 201(1 



Lokas et alj|2010l ; iKazantzidis et alj|201ll . and references therein 
for details). More recently accreted dwarfs (z < 1) may have a por- 
tion of their gas removed, with the remainder used in bursty star for- 
mation events as they pass close to the host throughout their orbit. 
Baryonic mass loss may be completed within a few million years, if 
the feedback is quite strong, or over hundreds of millions of years, 
as each feedback episode removes more mass and decreases the 
depth of the galaxy's potential well. 

If some dwarfs existed within groups before they were ac- 
creted by a primary halo, they could have undergone signif- 
icant additional mass loss. The existence of such a grouping 
befo re infall is under discussion for the Local Group dwarfs 
(seelD'Onghia & Lakdl2008l ; iLi & Helmill2008l ; iMetz et alj[200g| : 
IKlimentowski et alj I2OI0T " This means that the gas poor dSphs 
we currently observe may have undergone disruption for an 
even longer time period than expected. These interactions be- 
tween dwarfs can lead to a factor of two decrease in t> ma x 
and significant mass loss 1 Kravtsov, Gnedin, & Klypid 1 20041; 
Diemand. Kuhlen. & MadaJ 20071 ; iMadau. Diemand. & Kuhlerl 
20081 ; IKlimentowski et alj|2009i) . 

The inner density profile of the dwarf, the orbital pa- 
rameters, the inclusion of a baryonic disc in the host, and 
the mass ratio of dwarf to host galaxy all impact the severity 
of stripping that a satellite undergoes. Stronger inner den- 
sity cusps result in dwarfs that per sist even after undergoing 
the m ost e xtreme tidal stripping ( Colpi, Ma yer. & Governatol 
19991; IPefiarrubia. McConnachie & Nayarrol l2008at 



l20ld) . The inclusion of a disc component in the host can decrease 
the masses of satellites by an order of magnitude as well as the 
total number of subhaloes at redshift zero by a factor of two or 
three, depending on the ratio of host viria l mass to disc mass 
jPenarrubia elalfcoiol; iD'Onghia et"aill2010h . 

Much work has been done quantifying the effects of tidal 
stripping and heatin 



Tavlor & Babul 



Penarrubia et alj 



:ating (e.g. iGnedin, Hernquist. & Ostrikerl 1 1 99S 
200 li; iTaffoni et alj 120031; IRead et alj l2006bllc 



2001 ; I lattom et alj I2UU3I; IKead et alj !2UU6bUct 
2009l ; lD'Onghia et alJboiq) . Tidal heating is im- 



Penarrubia. Navarro. & McConnachie] 2008bt Penarrubia et alj 



portant for orbits with pericentres within 20 kpc, while disc shock- 
ing dominates for pericentres within 12 kpc. However, when a disc 
is present, orbits with pericentres within 20 kpc are completely de- 
stroyed, except those with the largest inner density profiles. The 
most massive satellites (M sa t /A/host > 0.1) sink quickly due to 
dynamical friction and merge with the host, whereas the least mas- 
sive satellites (M sa t/Mh oa t < 0.02) have no change in their orbital 
parameters but undergo significant mass loss. 

Both internal and environmental effects strongly influence the 
mass assembly history of a dwarf galaxy. This can result in dramat- 
ically different density profiles for galaxies of the same total mass. 
Due to the lack of certainty in how these haloes assembled, we test 
a variety of methods to remove the baryonic component of dSphs. 



3 METHODOLOGY 

3.1 Cosmological hydrodynamical simulation 

Realistic estimates of how a satellite's mass distribution changes 
due to the removal of baryons can be obtained by running cosmo- 
logical hydrodynamical simulations. To this extent, we ran cosmo- 
logical simulations of a dwarf galaxy with M v x t = 3 x 10 Mq 
with the same initial conditions twice, once without gas removal 
and the other time with the gas blown away just before 2 = 1. 

To perform the simulations, we used the adaptive 
mesh refinement A^-body+hydrodynamics code hydroART 
dKravtsov. Klypin. & Khokhlo\ll l997l ; f Kravtsovl 19991 ). The code is 
adaptive in both space and time, achieving higher resolution in re- 
gions of higher mass density. We used the following cosmological 
parameters: Q.m = 1 — S^a = 0.27, Sib = 0.04, and as = 0.8. 
The simulation zoomed in on a single halo inside a 10 /i _1 Mpc 
comoving box and achieved a maximum resolution of 60 pc and a 
DM particle mass m p ~ 10 5 Mq at redshift zero. 

The physical model used by the code includes many rele- 
vant processes, such as cooling down to 300 K from metals and 
molecules, a homogeneous UV background, metal advection, and 
stellar mass loss. Star formation occured stochastically in regions 
where the gas density reached a critical threshold, tisf > 1 cm -3 , 
with a constant efficiency, esf = 2.5 per cent. Each "star" par- 
ticle represented a simple stellar population with a Chabrier IMF 
and deposited thermal energy via SNe and stellar winds as well as 
radiation pressure from O and B stars into the surrounding ISM. 

This model is based on observations of star forming regions. 
More details on the implementation of the star formation and feed- 
back model can be found in Trujillo-Gomez, S. et al. (in prep). 
The isolated dwarf galaxies simulated using this new model are in 
very good agreement with observations. In particular, they have no 
bulges, slowly rising rotation curves, and constant or slowing in- 
creasing star formation histories. 

We created two isolated dwarf galaxies. The first run, the fidu- 
cial run, was performed as described above with no alterations. We 
also ran a second simulation, which we call the "no gas" run. It is 
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Table 1. Mass values for the hydroART simulated galaxies. Two isolated dwarfs were run, one fiducial model (left) and one where 
all gas within the galaxy was removed 200,000 years before z = 1 (right). The fiducial model increases its stellar mass and accretes 
more gas; however, the no gas dwarf retains only a hot diffuse gas component. The stellar mass in the no gas model is reduced from 
z = 1 to z = due to stellar mass loss, while the virial masses do not increase substantially from z = 1 to z = for either model. 



Model 


fiducial 






no gas 




M viT 


M gas (r < lOkpc) 


M* (r < 10 kpc) 


M vir 


M gas (r < 10 kpc) 


M* (r < 10 kpc) 


(M ) 


(M ) 


(M Q ) 


(M ) 


(M ) 


(M Q ) 


z = l 2.3 x 10 10 


6.7 x 10 s 


5.4 x 10 7 


2.2 x 10 10 


6.2 x 10 3 


5.4 x 10 7 


z = 3.0 x 10 10 


6.4 x 10 s 


3.0 x 10 8 


2.5 x 10 10 


6.2 x 10 3 


4.9 x 10 7 



the same as the fiducial run until just before z = 1, when we re- 
moved the gas in the model to emulate the gas stripping that occurs 
as a satellite falls into its host halo. To achieve this, we artificially 
increased the temperature of all the gas in the high-resolution re- 
gion to 10 6 K and halted cooling in those cells from z = 1 until 
z — 0. This removed the majority of the gas mass, with all bound 
gas existing as a hot diffuse medium. Only ~ 6000 Mq of gas re- 
mained in the galaxy with a total of ~ 10 6 Mq of gas within the 
virial radius of the halo. Since this gas could not cool, the galaxy 
was unable to form any more stars. 

These simulations do not include the effects of tidal stripping 
of DM and stars by the host potential, and as such, they provide 
an estimate of the minimum effect of baryonic physics on the DM 
density profiles. TableQjsummarises the properties of the two mod- 
els. In Section[4] we will use the results of these simulations to test 
our model for adiabatic expansion of the halo due to the baryonic 
loss. 



3.2 iV-body simulations 

We performed controlled collisionless iV-body simulations to 
quantify the gravitational effect of the removal of baryons from 
subhaloes and the enhanced tidal effects on subhaloes from the 
Milky way disc. In what follows, we describe the code and the ini- 
tial conditions we used. These simulations incorporated the undis- 
puted net results of baryonic processes. 



3.2.1 Code and initial conditions 

We used a direct summation iV-body code with the leapfrog time- 
stepping integration scheme, a constant time step At = 6.4 x 
10 yrs, and a Plummer softening of e — 20 pc. All resultant 
simulations retained an energy conservation of AE/E < 10 _J . 
The code was run in parallel by using both OMP (Open Multi- 
Processing) and MPI (Message Passing Interface) libraries. 

Initial conditions were set up by solving the equilibrium Jean's 
equation to obtain the rms velocity disp ersion « rms (r) as a func- 
tion of radius for a dwarf with a spherical iNavarro. Frenk. & White! 
dl996bl . fi 9971 . NFW profile hereafter) density distribution given by 



n{x) = ln(l + x 



AtrGr 2 fi(x max ) x(l + x) 2 ' 



r/r s 



(1 + x) 



(1) 

(2) 



Here « max is the maximum circular velocity of the dwarf and a; ma x 
is the distance at which the velocity peaks, where x max = 2.163. 

The distribution was truncated at the outer radius r ut = 16r s . 
This truncation radius is close to the expected virial radius for 
haloes with masses typical for large satellites of the Milky Way. 



The exact value of r out is unimportant because the satellite will be 
severely tidally stripped when it falls into the potential well of the 
MW. 

When setting initial conditions, we used multiple masses for 
particles. In the inner regions, the particle mass, mi, was small and 
nearly constant. The mass increased with distance as 



mi(r) = mi(0) 



1 + 



(3) 



where mi (0) is a normalisation constant. This prescription allowed 
us to resolve the central region of the dwarf with small particles and 
at the same time have an extended halo. Throughout the dwarf's 
evolution, particles of different mass moved in radius so that mas- 
sive particles occasionally came to the centre. We monitored the 
situation and found that the central ~ 1 kpc region was always 
dominated by very small particles with a very low contribution of 
more massive particles. As an example we quote results from a 
200,000 particle simulation after 10 dynamical times. Within the 
central 100 particles, the average particle mass increased by just 
2 per cent. Inside the radius of the first 1000 particles, the mass per 
particle increased by 6 per cent. 

We assumed an isotropic Gaussian distribution of velocities. 
Gaussian velocities produce small non-equilibrium effects for the 
NFW distribution. We allowed the system to evolve in isolation 
for 10 dynamical times to settle and come into equilibrium. Non- 
equilibrium effects were small at any moment. For example, af- 
ter 10 dynamical times, the maximum circular velocity was within 
2 per cent of initial value. The largest effect was a decline in the 
density, which occurred at the very centre. The density in the cen- 
tral ~ 100 pc declined by about 30 per cent. We used the evolved 
configuration as the reference when we examined the results of 
other simulations. Overall, these changes were at a low enough 
level for the density and circular velocities to retain their basic 
NFW shape and cusp. 

Our simulations consisted of 200,000 particles of various 
mass. If we had used a constant mass per particle, we would have 
needed N e g = 1.3 X 10 6 particles to resolve the central region with 
the same small particles. 

The gravitational force of the MW acting on satellites was 
modelled using frozen external potentials, which roughly approxi- 
mated the mass distribution in our MW galaxy. 



3.2.2 Scaling models 

When scaling TV-body simulations to physical units, we chose the 
radius (ro) and velocity (vo) to be our free parameters, which left all 
other quantities, such as time-scale (to) and mass (mo), dependent 
on the choice of radius and velocity scales. 

Our simulated dwarfs were created to test the effects of baryon 
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removal and tidal stripping for large satellites. The initial dwarf has 
the following properties: r s = 4.0 kpc and w max = 62 kms -1 . 
This maximum circular velocity is somewhat high for dSphs. How- 
ever, these parame ters were chosen to mimic the high end of the 
"massive failures" jBovlan-Kolchin et all20l"3) before infall. Sim- 
ilar to the results of iBovlan-Kolchin et alj l l2012h . the entire initial 
profile of our dwarf is inconsistent at 2a with all the dSphs in the 
Local Group. This initial setup allows us to test if an extreme case 
for a dwarf satellite is able to reproduce observations by including 
baryonic effects. 

Our host galaxy was matched to the Aquarius E galaxy 
virial dark matter mass with an added baryonic component. 
The host galaxy is denned as: rdi sc = 3 kpc, Mdisc = 
6.05 x 10 10 M , r s = 25.12 kpc and M vi[ = 1.39 x 
10 12 M (i; circ = 179 kms" 1 ). These are reasonable val- 
ues for the MW halo and disc mass, which are estimated to 



be Mtot = 0.7-3 x 10 12 iwilkinson & Evans! 1 19991; 



Sakamoto, Chiba 


& Beersl 120031: ISmith et ail 120071: iLi & White! 


20081; Xueetal. 


20081; Kallivavalil et al 20091: Guo et alj 20ld: 


Watkins, Evans & An! 2010b and M diar ~ 5.5 x 10 10 


jGerhard 2002; Flvnn et al.l 2006). respectively. 



We calculated the circular velocity of the initial perturber at 
8 kpc, to provide the reader with a comparison of the halo and the 
MW. Using the mass within a given radius equation for the expo- 
nential disc and the NFW profile, we found that the baryon values 
are M hB , r (r < 8 kpc) = 4.5 x 10 10 M or a « circ = 156 km s" 1 
and the DM values are M DM (r < 8 kpc) = 4.9 x 10 10 M 
or a v c i IC = 162 km s _1 . After summing the circular velocities 
in quadrature we obtained i; c i rc = 224 kms -1 . This is close 
to the circular velocity of th e solar neighbour hood v c irc(-Ro) = 
218 ± 6 km s" 1 , as found bv lBovv et ait < l2012l) . 

3.2.3 Setup of tidal stripping simulations 

The simulations were run in the presence of an external MW-sized 
perturber after running in isolation. The no mass removal case (NR) 
was run with a pure NFW perturber, while all other cases were 
run with a NFW plus spherical exponential "disc" perturber. This 
setup allows the NR case to remain as the pure DM example, while 
the mass removal cases show the effects of including baryons (by 
adding a disc) in the simulation (see Table [2] for naming conven- 
tions). 

All satellites were run in circular orbits at r or bit = 50 kpc, 
70 kpc, 100 kpc, and 150 kpc to examine the most extreme effects 
of tidal stripping from different distances. Elliptical orbits were also 
run to give more realistic results of tidal stripping. We created orbits 
with varying eccentricities by changing apocentre and pericentre 
distances from 1:2 to 1:5 to compare the effects of different sub- 
halo orbital evolutionary histories. These or bits are less e l liptica l 
than the typical orbital ratio of 1 :6 derived bv lGhigna et~aT1 dl998h ; 
however, they are we ll within the 90 per cent interval found by 
iDiemand et alj J20071) . which extend s from ^1:28 to ^2:3 orbits. 

In the Via Lactea simulations jDiemand et alj|2007h . almost 
all galaxies within the "earliest forming substructure" sample had 
orbits that fell within 30 kp c of the host. In the SPH simulations 
run bv lZolotov et alj d2012l) . over half of all the subhaloes had or- 
bits with pericentres below 30 kpc (private communication). In our 
truncation radius estimates (Section I5.lt , we find that a satellite 
in the presence of a MW with a baryonic disc component should 
have interesting tidal evolution if it exists in a orbit with pericentre 
~ 40 kpc. This occurred when the truncation radius, r t , had values 
from 2 to 4 times the scale radius, r s . We ran four elliptical orbits. 



Table 2. Different dwarf galaxy cases. The initial and final maximum cir- 
cular velocities during their isolated runs are listed. 



Label 


Description 




^max,f 


"max.f 
^max i 


NR 


no baryonic mass removal 


62.1 


64.2 


1.03 


PR 


pre-halo formation removal 


64.2 


57.6 


0.90 


ER 


exponential removal 


64.2 


51.2 


0.80 


IR 


instantaneous removal 


64.2 


51.0 


0.79 



Table 3. Orbital parameters for each satellite galaxy. Included are their ini- 
tial maximum circular velocities and final maximum circular velocities after 
5 Gyrs. Labels that end in "c" denote a circular orbit and those that end in 
"e" denote an elliptical orbit. 



Label 


Tlf^cr^finti rMi 


^max,i 




"max.f 
^max i 


NR150c 


NR, 150 kpc orbit 


64.2 


60.7 


0.95 


NRlOOc 


NR, 100 kpc orbit 


64.2 


57.2 


0.89 


NR70c 


NR, 70 kpc orbit 


64.2 


51.6 


0.80 


NR50c 


NR, 50 kpc orbit 


64.2 


43.0 


0.67 


NRl:2e 


NR, 150 to 70 kpc orbit 


64.2 


56.5 


0.88 


NRl:3e 


NR, 150 to 50 kpc orbit 


64.2 


53.8 


0.84 


NRl:5e 


NR, 150 to 30 kpc orbit 


64.2 


51.2 


0.80 


NRl:4e 


NR, 120 to 30 kpc orbit 


64.2 


46.9 


0.73 


PR 150c 


PR, 150 kpc orbit 


57.6 


53.6 


0.93 


PR 100c 


PR, 100 kpc orbit 


57.6 


49.5 


0.86 


PR70c 


PR, 70 kpc orbit 


57.6 


42.1 


0.73 


PR50c 


PR, 50 kpc orbit 


57.6 


28.1 


0.49 


PRl:2e 


PR, 150 to 70 kpc orbit 


57.6 


48.5 


0.84 


PRl:3e 


PR, 150 to 50 kpc orbit 


57.6 


44.7 


0.78 


PRl:5e 


PR, 150 to 30 kpc orbit 


57.6 


40.4 


0.70 


PRl:4e 


PR, 120 to 30 kpc orbit 


57.6 


34.2 


0.59 


ER150c 


ER, 150 kpc orbit 


51.2 


46.0 


0.90 


ERIOOc 


ER, 100 kpc orbit 


51.2 


40.5 


0.79 


ER70c 


ER, 70 kpc orbit 


51.2 


30.0 


0.59 


ER50c 


ER, 50 kpc orbit 


51.2 


10.7 


0.21 


ERl:2e 


ER, 150 to 70 kpc orbit 


51.2 


39.6 


0.77 


ERl:3e 


ER, 150 to 50 kpc orbit 


51.2 


34.8 


0.68 


ERl:5e 


ER, 150 to 30 kpc orbit 


51.2 


29.9 


0.58 


ERl:4e 


ER, 120 to 30 kpc orbit 


51.2 


22.5 


0.44 


IR150c 


IR, 150 kpc orbit 


51.0 


45.8 


0.90 


IRlOOc 


IR, 100 kpc orbit 


51.0 


39.9 


0.78 


IR70c 


IR, 70 kpc orbit 


51.0 


29.0 


0.57 


IR50c 


IR, 50 kpc orbit 


51.0 


10.6 


0.21 


IRl:2e 


IR, 150 to 70 kpc orbit 


51.0 


39.0 


0.76 


IRl:3e 


IR, 150 to 50 kpc orbit 


51.0 


34.0 


0.67 


IRl:5e 


IR, 150 to 30 kpc orbit 


51.0 


29.2 


0.57 


IRl:4e 


IR, 120 to 30 kpc orbit 


51.0 


21.4 


0.42 



Three of the four orbits had the same apocentre and are listed here: 
l:2e (150-70 kpc), l:3e (150-50 kpc), and l:5e (150-30 kpc). The 
fourth orbit, l:4e (120-30 kpc), had a different ellipticity but ar- 
rived at the same 30 kpc pericentre as the l:5e orbit. A description 
of each dwarf galaxy's orbits is provided in Table|3] 

3.3 Adiabatic expansion 

The removal of baryons within a galaxy results in smaller mass 
inside a given radius r. In turn, this reduces the force of grav- 
ity and causes the galaxy to expand. This is a known and qual- 
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itatively well understood physical process jZeldovichetalJll98d : 
lBlumenthaletal.lll98d iGnedin etal]|2004ir It is also observed 
in many cosmological simulations that include hydrodynam- 
ics and star formation (IGustafss on. Fairbairn, & Sommer-Larsen 



20061: IColm Valenzuela. & Klypinl l200d : iTissera et alj 1201a 
Gnedin et alJuOllh. He re we use the simplest prescription of 
Blumenthal et al .1 a 1986T) to make an analytical model for the ex- 
pansion of a galaxy when it loses mass. It assumes circular orbits 
for dark matter particles and uses the conservation of angular mo- 
mentum to derive a relation between the initial mass, Mi(n), and 
radius, n, and the final mass, Mf(rf), and radius, rj. Note that 
the prescription does not use adiabatic invariants; it uses an in- 
tegral of motion - the angular momentum. Thus its limitation is 
not that the process must be slow (adiabatic), it is the assumption 
of circular motions. An extension of the prescription, provided by 
IGnedin et alj 1 2004 ), makes a correction due to non-circular mo- 



Gnedin et alj J2004) shows, in the regime of small com- 



tions. As 

pression (or expansion) the difference between the two approxima- 
tions is small. 

It is convenient to rewrite the iBlumenthal et alj Jl986l) pre- 
scription using the circular velocity, V c irc = y GM(r) /r of a 
mass shell, where M(r) is the mass inside radius r. If a(ri) is the 
fraction of mass that remains inside n after the removal of baryons, 
then the final radius, rt, and circular velocity, Vi , are given by 



Vf{rf) — a(r-i) Vj(n) = a(ar{) Vi(art). 



(4) 
(5) 



Note that the factor of a in equation (|5j enters linearly in both 
the radius, ri, and circular velocity, Vi. As a result, the effect of 
the expansion is rather large. This should be contrasted with the 
case that lacks expansion when mass M(r) decreases by a factor 
of a but the radius stays the same: r-f = r^. In this case, the fi- 
nal circular velocity, Vf, declines only by a factor of ^fa, which 
is a small effect. For the baryon removal case, where /b ~ 0.17 
dKomatsu et~ai]l201ll) . we find a ~ 10 per cent change, where it is 
expected that a ~ 1 — /b. 



3.4 A tidal stripping model 

Simulations of the effects of tidal stripping can be very demand- 
ing and take substantial computing resources. It is helpful to have a 
simple model that can reproduce some of the simulation results. We 
found that the following model produces quite accurate results. The 
motivation is simple; after the removal of some fraction of mass 
due to tidal forces, the halo expands and some particles may be- 
come unbound giving rise to another round of mass loss. Note that 
this happens because dark matter particles move on elliptical or- 
bits, which are typically very elongated. Thus the removal of mass 
from the peripheral regions, where the tidal force is strong, leads to 
a lack of mass in the inner regions over time. Therefore, even the 
central region of a halo may be affected by the tidal stripping. 

A realisation of the NFW distribution was made using the pre- 
scription described in Section 13.2.1 1 For this model we increased 
the number of particles to 500,000, which provided many low mass 
particles in the central region of the dwarf. This allowed us to make 
detailed calculations of tidal stripping without numerical errors. 

After the dwarf was created, we estimated the energy of each 
particle and removed all unbound particles. When estimating the 
potential energy of a particle, we included the potential associated 
with a spherically symmetric tidal force. Only particles inside the 



tidal radius were used for calculations of the gravitational potential 
of the dark matter itself. 

The initial tidal radius, r t , was parametrised by its ratio to 
the characteristic radius, r s . This set the normalisation for the tidal 
potential, which did not change when we proceeded with iterations 
of tidal stripping. 

Once the unbound particles were removed, the halo was adia- 
batically expanded following equations Q and {5). Following this 
method, we found the mass of bound particles within each radius 
and compared it with the mass at the previous iteration. The ratio 
of these masses gave the a parameter, which was then used for the 
adiabatic expansion; the coordinates were increased and the veloc- 
ities were decreased according to equations l(4) and (|5). We also 
found the new tidal radius and proceeded with the next round of 
removal of unbound particles. 

Just like in our iV-body simulations, we found that tidal strip- 
ping continuously removed mass from the halo. In other words, 
there was no convergence, and some fraction of mass was removed 
even after many iterations. However, the rate of the mass removal 
dramatically depended on the initial tidal radius. We found very lit- 
tle effect when the initial tidal radius exceeded r t = 4r s . At smaller 
tidal radii, the rate of stripping increased. By comparing the simple 
model to the TV-body simulations, we found that between 10 and 
30 iterations gave results very similar to those observed in actual 
simulations. 



4 EFFECTS OF BARYON REMOVAL 

4.1 Effects of baryon removal in hydrodynamical 

simulations: testing the adiabatic expansion model 

In this section, we examine the effects of infall and gas removal on 
the dark matter mass distribution of a dwarf galaxy. We modelled 
this process in detail by performing two high resolution cosmolog- 
ical simulations of the formation and evolution of a dwarf galaxy, 
including gas, star formation, and feedback. 

The fiducial run was done with no removal of gas - the system 
evolved without any intervention - and mimics the evolution of a 
dlrr galaxy. As the left panel of Fig.Q]demonstrates, the simulation 
reproduces some generic properties of dlrrs: almost no bulge, gen- 
tly rising rotation curve and a large mass in gas as compared with 
the stellar mass. The model does have an issue related with the 
strong stellar feedback used for the simulation; the fraction of mass 
in baryons is a factor of two too low. In this simulation, there are 
~ 0.5/b the central region, r < 5 kpc, while obse rvations point t o 
a larger value, closer to the cosmological fraction (lOh et al JI201 lb . 
In the second run, the gas was removed from the halo just before 
z — 1 and evolved to z = 0, without allowing the gas component 
to cool. 

Note that the total circular velocity profile of the fiducial run 
does not change appreciably from z = 1 to z = 0, without inter- 
vention. There was little evolution in total mass in the outer regions 
of the dwarf; in contrast, within the inner 500 pc, the mass de- 
creased slightly. The lack of evolution in the fiducial isolated dwarf 
over this time period, with the exception of continued star forma- 
tion, makes this particular halo a good test of the adiabatic expan- 
sion model against a simulation with gas removal. 

Fig. compares the fiducial and no gas removal models. Re- 
moval of gas clearly resulted in the expansion of the halo, with 
the circular velocity declining at every radius. The fraction of mass 
removed at z — 1 is small; within the central 5 kpc region, the 
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Figure 1. Results from the cosmological hydrodynamical simulations displaying the effects of gas removal from an isolated dwarf galaxy. Left: The circular 
velocity profile of the fiducial model at Z = 1 is shown as the solid line. The dotted (red) curve shows the gas removal simulation results at z = 0. The 
short-dashed (blue) curve is for the adiabatic expansion model. The bottom curves represent the circular velocity profiles of the gas (dot-dashed curve) and 
stars (long-dashed curve) in the simulations at z = 1 before the gas removal. Right: Here, the ratio of total mass within a given radius at z = over the total 
mass within a given radius before gas removal (at z = 1) is shown. The dotted (red) curve shows the gas removal simulation results, and the short-dashed 
(blue) curve gives the predictions of the adiabatic expansion model. The model gives 10 per cent accuracy for radii r > 500 pc and underpredicts the effect at 
smaller distances. 



fraction of mass removed is Mgas/Mtot = 0.09. It is even smaller 
(0.02) inside the central 300 pc region. However, the final expan- 
sion was substantial. For example, the mass inside ~ 500 pc de- 
creased by ~ 40 per cent. 

The results of the cosmological hydro simulations, in agree- 
ment with results presented later from iV-body simulations, show 
that the expansion of the halo takes place when a fraction of mass is 
removed from the central region. We approximated removal of the 
baryonic component of the hydro simulations in two ways. In our 
first trial, we removed the fraction of gas mass wit hin a given radius 
from t he dwarf's total mass and implemented the lBlumenthal et al] 
expansion. This fraction varies with radius and produced 
weak expansion at small radii, as there was little gas mass within 
central ~ 1 kpc. In our second trial, we removed a constant mass 
fraction of 10 per cent at all radii, which corresponded to the aver- 
age baryonic mass lost within the entire dwarf. Our second method, 
using constant removal, produced results in better agreement with 
the no gas hydro sim ulation. 

In this case, the lBlumenthal et"al prescription (equa- 

tions [4] and [5} underestimated the expansion in the very centre but 
gave reasonably accurate results for radii larger tha n ~ 0.5 kpc. 
Const ant mass removal is also more in line with the iGnedin et alj 
d2004l) model. They found that expansion corresponds to more 
mass removed within a given radius than the expected value, after 
accounting for non-circular motions. The analytical model has the 
advantage that it can be used for different parameters and different 
configurations, so we use it to investigate the effects for different 
amounts of mass removal. 

We consider a simple but realistic model in which baryonic ef- 
fects cause the initial halo to retain a constant fraction of its density 
at every radius: a — ao = const. In this case, adiabatic expansion 
of the halo preserves the shape of the density distribution, but with 
modifications to its parameters. Any halo profile can be written in 
the form 



p(r) = po* — 
So 



where ro is a scale radius and po defines overall normalisation of 
the profile. The dimensionless function ty(x) describes the density 
profile. For example, for the NFW profile (equations [JJ and [2] one 
can chose ro = r s and ^ = l/x(l + x)' 2 . When the halo loses 
1 — ao of its mass, it expands and changes the parameters ro and 
po, but the function is preserved. Using equations (|4j-((6) and 
labelling variables with the subscripts i for initial and f for final 
parameters, we get 



po,f = Qo P0,i 



n>,f = 



ro,i 



(7) 



Note the large power of ao in the equation for density normalisa- 
tion. It means that the overall normalisation for the density declines 
very substantially. If we assume that most of the baryons were 
removed when a din galaxy transforms into a dSph, then ao ~ 
1 — /b ~ 0.8, which leads to a drop in density of po,f / po,i ~ 0.4. 
Next, we pose the question of how much the density changes at a 
given physical radius after baryon removal. The answer depends on 
the particular form of the density profile: 



Mr) 
pi(r) 



4 

= «0 



V r o,f / 



(8) 



If the function <5 is the power law oc x , then 

Mf) 4-/3 



(9) 



(6) 



For the NFW profile, the central cusp has a slope of j3 — 1 and the 
density declines by large factor, pi/pi — a'o ~ 0.5, again assuming 
that ao ~ 0.8. The change in the outer region is much smaller as 
ft — 3, so pi/ pi — ao ~ 0.8. It is interesting to note that a density 
profile with a central core (/3 = 0) experiences the largest decline 
in density, pij p\ = a%. 
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Figure 2. Profiles of dwarf galaxies at their final moment in isolated evolution. In all plots, the short-dashed (black) line shows the no mass removal case 
(NR), the dash-dotted (red) line shows the mass removal before the halo formed (PR), the solid (blue) line shows the exponential removal case (ER) and the 
long-dashed (grey) line shows the instantaneous removal case (IR). Top Left: The circular velocity profiles of the four dwarf galaxies. An additional line, dotted 
(purple), displays the adiabatic expansion model that the NR halo would incur if it lost 20 per cent of its mass. It lies on top of the ER and IR profiles, indicating 
they underwent adiabatic expansion. Bottom Left: The density profiles of all four dwarfs over the same range as above. Top Right: The density profiles of all 
four dwarfs, zoomed into the central 50 pc to 1 kpc region. The initial state of the dwarf, before any isolated evolution or mass removal, is included as the 
dotted (purple) line. The NR halo begins to deviate from this profile at 500 pc. Bottom Right: The fractional deviation in density of the mass removal cases 
divided by the NR profile. The profiles in the right two panels represent the mean of 21 profiles in order to reduce shot noise. The errors displayed are the 
associated errors in the mean. 



4.2 Baryon removal in iV-body simulations 

The correct ratio of baryons to dark matter is important for an 
accurate representation of the structure of dwarf galaxies. Dwarf 
spheroidals are known for their high mass to light ratios, producing 
distributions that do not follow the cosmic mass fractions. In order 
to mimic this lack of baryons inside of Sphs, we removed some 
fraction of mass from our dwarfs. The mass removed should reflect 
the amount of baryonic matter lost when accounting for the lack 
of gas in dSphs. This left the DM mass of the galaxy unchanged. 
As the details of removal process are still unknown, we simply re- 
moved some fraction of the mass from our test dwarf in isolation 
to create a comparison between dwarfs from pure DM simulations 
with NFW profiles and those with baryons included. 

To test a few cases, we have run four dwarfs in isolation under 
different conditions (see Table|2]for naming conventions). Our test 
dwarf (no mass removal - NR) was run for 1 Gyr from a equilib- 
rium setup and is used as a benchmark for our comparisons to other 
dwarfs. 



The three remaining dwarfs were modified using techniques 
that test a variety of ways a dwarf galaxy may have its baryonic 
component removed. On one extreme, we present the case where 
the majority of baryons have never entered the dark matter halo. 
Here, the dwarf will always be in equilibrium with the correct mass 
fraction, and there is no need to expel more baryons. This creates 
a steep inner density profile making mass removal via tidal strip- 
ping difficult. We explore this scenario with a dwarf that undergoes 
pre-halo formation mass removal (PR). It was evolved in isolation 
with an initial mass at 80 per cent of the no mass removal case and 
remained in equilibrium. 

However, it is more likely that the dwarf once contained a 
baryonic component that was either expelled via SNe feedback, 
galactic winds, and radiation pressure or removed by ram pressure 
and tidal stripping. One dwarf (instantaneous removal - IR) began 
as the no removal case and then had 20 per cent of its mass instan- 
taneously removed before evolving in isolation. 

Another dwarf (exponential removal - ER) explored slow 
mass removal throughout its evolution in isolation. The mass be- 
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gan at the no removal level and then was reduced via an expo- 
nential function of the form: mi(t) = mi_fi{B + Ae~ t ^ t °). Each 
particle began with a mass of mj,o and currently has a mass of 
mi(t). We have implemented to ~ 200 Myrs with B = 0.8 and 
A — 0.2. This equation allowed 20 per cent of the mass of each par- 
ticle to be removed over an exponential time-scale, to. The larger 
the value of the time-scale, the slower mass is removed. The re- 
moval time-scale should be larger than the dynamical time of the 
central region. Using tdyn = R/vdrc we found t<x yn ~ 23 Myrs 
for R — 1 kpc and t; c i rc = 43 km s~ and tdyn ~ 16 Myrs for 
R = 0.5 kpc and v c i IC — 30 km s _1 . It is apparent that the exact 
value depends on the definition of the inner region. Our choice of 
time-scale, to ~ 200 Myrs, is an order of magnitude larger than 
these calculations. 

In both the IR and ER cases, the galaxy initially had the same 
mass as the no mass removal case. When the baryons are removed, 
the halo is expected to adjust to the mass loss according to adia- 
batic expansion. The PR case examines how a dwarf of the same 
mass without expansion evolves. This tests the differences between 
haloes that undergo baryon removal with and without expansion. 



4.2.1 Effects of fast and slow loss of baryons 

The change in maximum circular velocity, u max , throughout each 
dwarf galaxy's isolated evolution is given in Table [2] For the no 
mass removal case (NR), the initial n max in the table is from the 
equilibrium conditions given by the code before any isolated evo- 
lution, and the final u max is after isolated evolution. This shows 
how much the initial profile we provide relaxed after 10 dynamical 
times. For all other cases, the initial w max listed is the final state of 
the NR case to show the deviation from the preferred equilibrium 
state of the dwarf galaxy. 

Fig. [2] shows the density and circular velocity profiles of the 
dwarfs at their final moment of isolated evolution. We see that the 
removal of baryons after the halo has formed, via instantaneous 
mass removal (IR) or exponential mass removal (ER), dramatically 
reduced the circular velocity profile. Mass removal caused the halo 
to adiabatically expand, which reduced « max and increased r max , 
the radius where the velocity peaks. This can be seen in the differ- 
ence between the NR halo and the ER and IR haloes. Our adiabatic 
expansion model, explained in Section [33l is also plotted as a cir- 
cular velocity profile in Fig. [2] It began with the initial mass profile 
of the NR halo and used an a = 0.8 to reduce the v max and in- 
crease r max . We find that our simple adiabatic expansion model re- 
produces the final circular velocity profiles of the ER and IR haloes 
very well. 

In contrast, the removal of baryons before the halo has formed, 
as done in the PR case, caused u max to decrease without an increase 
in r max . This is because the halo did not undergo adiabatic expan- 
sion, as there is no baryonic removal from the halo, but instead the 
baryons were prevented from ever falling into the halo. This cre- 
ated a higher central density profile with a more pronounced cusp 
than in the adiabatic expansion cases. 

The density profiles of the four dwarfs are also shown, with 
the largest difference between the simulations located in the central 
1 kpc. To better examine the central density of the dwarfs, the right 
two plots in Fig. [2] zoom into the inner 50 pc to 1 kpc region. The 
upper right plot in Fig.|2]shows inner density profiles and includes 
the original NFW halo without mass removal and before relaxation. 
The NR density profile has relaxed for 1 Gyr, and differences be- 
tween the initial NFW halo and the NR halo are only seen within 
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Figure 3. Dependence of the satellite truncation radius (in units of the scale 
radius) on the radial distance from the centre of its host halo. Satellites of 
varying masses (M 8a t = 10 10 , 3 X 10 10 , and 10 11 Mq) are shown from 
top to bottom and in black, red, and blue, respectively. Dashed lines show 
the truncation radii for satellites orbiting in the presence of a pure NFW halo 
(M v ; r = 1.4 X 10 12 Mq) as contrasted with solid lines of the same colour 
showing a primary halo that also contains baryons (iWdiso = 6 X 10 10 Mq) 
that coalesced at the centre of the halo. The horizontal dotted line denotes 
rt = 2.2r s . The internal structure of subhaloes can be dramatically affected 
when their orbits include pericentric passages with truncation radii below 
2.2r s . 

the inner 500 pc. We compare all isolated dwarf cases to the NR 
case, as it represents the preferred equilibrium state of the dwarf. 

The lower right panel of Fig. |2]shows the density of the three 
mass removal cases divided by the NR case. Here, we examine 
how the density ratio differs for different mass removal cases and 
how this compares to a DM only dwarf. Both the IR and ER cases 
evolved to have similar fractional density profiles with only half 
the original halo's density in the central region. The halo that never 
allowed baryonic infall (PR) retained 0.8 of the NR density pro- 
file. All profiles remained fairly flat without a large central density 
deficit. A decrease would have indicated the creation of a core from 
the traditional cuspy interior of the NFW profile. However, we do 
not see this effect occurring as a result of mass removal with or 
without adiabatic expansion. 

The PR case will examine how tidal stripping impacts a satel- 
lite that does not experience expansion. The IR and ER cases ex- 
plore the combined effects of expansion and tidal stripping for a 
massive dwarf. We expect that the IR and ER cases will be more 
susceptible to tidal stripping than the PR case as they have experi- 
enced an increase in r max . 



5 EFFECTS OF TIDAL STRIPPING 

5.1 Truncation radius: analytical estimates 

In order to examine how the concentration of baryons in the inner 
part of the MW halo impacts the subhalo, we begin by presenting 
analytical estimates of the subhalo's tidal radius, depending on its 
mass, the main halo mass, its distance from the main halo, and the 
inclusion or absence of baryons in the centre of the primary. One 
needs to be aware that the fate of the dwarf strongly depends on 
the tidal radius, r t , to which it is stripped. We later show that if 
r\ > 4r s , there is no significant evolution of i) max - However, tidal 
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stripping can even alter the very centre of the dwarf when r t < 
2r s ~ r max . So there is a transition at r t ~ 2r s with a small effect 
for larger rt and drastic stripping at smaller tidal radii. 

If the subhalo has a mass profile of m(r), which orbits a halo 
with mass profile M(R) at a radius R, then the tidal radius for the 
subhalo, r t idc, is defined as the radius of the subhalo at which the 
gravitational force on a test particle from the subhalo equals the 
difference of the gravitational force (tid al force) from the h alo on 
the subhalo and the test particle (see e.g. lKlypinetalJll999j) , 



G 



m(rtidc) 
r 2 

'tide 

m(r tide ) 



_d_ 

dR 
M(R) 
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GM(R) 



R? 



nide 



R dM 
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(10) 



(11) 



The subhalo can also be stripped to an even smaller radius defined 
by resonances between the gravitational force from the subhalo and 
the tidal force from the halo. This resonant radius can be obtained 
by solving the following equation: 



m(r rcs ) = M{R) 
r 3 ~ R s 



(12) 



The smaller of rtide and r re s is taken as the truncation radius. This 
is denoted by r t and defined as the radius to which the subhalo is 
stripped. 

For haloes following the NFW density profile, the mass profile 
m(r) is given by 



m(r) = Mvh 



M( c vir) 



r 

X = — , 



(13) 



where r s is the scale radius of the halo and is related to the virial 
radius of the halo through the concentration parameter, c vlT = 
r v ir/»" s , and the function fi(x) is given by equation (|2}. We assume 
that the subhalo mass profile follows the above distribution. 

The main halo includes a dark matter component that follows 
the NFW profile and a baryonic component that is concentrated at 
the centre of the halo and is well inside the orbit of the satellite. 
Therefore, the total mass profile of the halo is given by 



M(H) = M bar + M vir (l-/ b ) 



Ai(Cvir' 



(14) 



where x = R/R v [ r and CVir is the concentration of the halo. The 
parameter / b is the fraction of the virial mass that is in the baryonic 
component: / b = fibar/fimatter. Here, fi bar and f) ma tter are the 
contributions of the baryons and the total matter as compared to 
the critical density of the Universe, respectively. We assume that the 
baryonic mass of the halo is M bar = 6 x 10 10 M^p , consistent with 
the baryonic mass estimates for the MW (e.g.. iDehnen & Binnevl 
ll998l : lKlvpin. Zhao. & Somervillell 20021) . 

To determine the truncation radius, r t , we substitute the mass 
profiles for the subhalo (equation 113) and the main halo (equa- 
tion Q3} into the tidal and resonant radii equations i l I 1 1 and j 1 2b . 
A range of subhalo masses and virial host masses are used to ex- 
amine how the mass of both the primary and secondary impacts 
tidal stripping. 

The results of these calculations are in Fig. [5] which depicts 
the dependence of the satellite truncation radius (in units of the 
scale radius) on the radial distance from the centre of its host halo. 
Decreasing values of r t /r s indicate more tidal stripping because, 
as the tidal radius approaches the scale radius of the satellite, an in- 
creasing amount of mass is stripped. The maximum circular veloc- 
ity occurs at a scale radius r s = 2.2; so when r t /r s ~ 2, internal 
structure is affected. Satellites of varying masses are shown with 




Figure 4. Evolution of the circular velocity profiles for two mass removal 
satellites after experiencing 50 kpc circular orbits. The exponential mass 
removal (ER) satellite is shown in solid lines. The upper and lower black 
lines indicate its initial and final states, respectively, while the grey lines in 
between are each separated by one Gyr. After 5 Gyrs, the satellite's Umax 
has been stripped from 50 to 10 km s _1 . The initial and final state of the 
instantaneous mass removal (IR) case is shown as dashed (black) lines for 
comparison. The black squares are observational data points of the half- 
light velocities for the classical Local Group dSphs. 



different colours. Smaller haloes collapse earlier, increasing their 
concentration. This results in a larger density at r s , causing the less 
massive satellites to be more resilient against tidal stripping. 

Fig. [5] includes results from a NFW main halo (M V1I = 
1.4 x 10 12 Mq) with and without a fraction of the baryons 
(Mdisc = 6 x 10 10 Mq) coalesced at the centre. The pure-NFW 
halo mimics cosmological simulations, where baryons are assumed 
to trace the dark matter in an unbiased way. The impact of including 
a disc component is clearly seen for orbits with pericentres within 
~ 30 kpc. At this radius, each satellite's tidal radius decreases by 
~ 14 per cent, which is a significant effect. For orbits with peri- 
centres within ~ 20 kpc, the satellites will be destroyed under 
the influence of a MW with or without a disc component because 
r t /r s < 1. In contrast, if the pericentre is greater than 100 kpc, 
then r t /r s > 4 regardless of whether a disc component is included 
or not, which produces weak tidal stripping. 

For the regime of orbits with 30 kpc < r < 100 kpc, the 
inclusion of a disc has a more subtle but important effect. From 
30 kpc < r < 60 kpc the satellites have obviously strong tidal 
stripping if the MW disc is included. Orbits with 60 kpc < r < 
100 kpc may also undergo increased tidal stripping and should be 
examined further. We predict that the inclusion of a disc for orbits 
with these pericentres results in an important regime of mass loss 
and circular velocity decrease. As such, we have created satellites 
on elliptical orbits with pericentres between 30 kpc and 70 kpc, to 
investigate this region. 



5.2 Effects of tidal striping in A-body simulations 

Each of our three dwarf galaxy mass removal cases had 20 per cent 
of its mass removed, but they produced different maximum circu- 
lar velocity values and central density profiles. Instantaneous mass 
removal (IR) and exponential mass removal (ER) scenarios both 
resulted in a 20 per cent reduction in « max due to adiabatic expan- 
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Fi gure 5. Comparison of our ci rcular orbits to observations of the MW satellites (black squares) and the Aquarius E halo's "massive failures" from the analysis 
of Bovlan-Kolchi n~t alj 1201 2h (thick grey lines). Left: The circular velocity profiles of the no mass removal case (NR). This profile does not agree with the 
MW dSph population and was created to mimic the Aquarius cosmological simulations. Right: The circular velocity profiles of the exponential mass removal 
case (ER). These profiles span a much larger range of parameter space and are in agreement with observed MW dwarfs due to their inclusion of baryonic 
effects. Both panels show the initial isolated profile as a solid black line. The circular orbits after 5 Gyrs of evolution are shown from top to bottom as: 150 kpc 
(dotted), 100 kpc (short-dashed), 70 kpc (dash-dotted), and 50 kpc (long-dashed). 



sion. The pre-halo formation mass removal scenario (PR) did not 
undergo adiabatic expansion; therefore, its u max only decreased by 
10 per cent in isolation. See Section l4~2l for the mass removal de- 
scription and results. 

After each dwarf evolved in isolation, they were put into a 
variety of circular and elliptical orbits around a MW-like host as 
described in Section 13.2.31 and listed in Table [5] Comparing the 
orbital evolution of these three haloes will provide insight into how 
pericentre passage, number of orbits completed, and central density 
of the satellite impacts the final maximum circular velocity values. 



5.2.1 Circular orbits 

Four different dwarfs were placed in circular orbits around the MW 
perturber. The no mass removal dwarf (NR), which had evolved in 
isolation, was placed around a pure-NFW perturber to mimic the 
results of DM only simulations. The three mass removal dwarfs - 
two that showed adiabatic expansion (IR and ER) and one without 
(PR) - were placed in circular orbits around a MW perturber with 
a NFW and spheroidal exponential "disc" component. All satellites 
were simulated for 5 Gyrs with initial and final maximum circular 
velocities listed in Table[3] 

We investigated two regimes (i) at r ^ 100 kpc where the 
mass loss was expected to be low based on our truncation radius 
analysis (r t /r s > 4) and (ii) at 30 kpc < r < 100 kpc where the 
truncation radius ratio had values between 2 < r t /r s < 5. This 
second regime is a sensitive area where a large amount of mass 
may be lost if the tidal radius becomes small enough. For pericen- 
tre passages within 30 kpc, r t /r s < 2 resulting in extreme tidal 
stripping for the satellite regardless of its mass or central density. 

To investigate the first regime, two different circular orbits 
were created. One was created at 150 kpc and another was cre- 
ated at 100 kpc. The truncation radius to scale radius ratio (r t /r s ) 
is large at this orbital distance (for NR with a pure NFW MW: 
r t /r s = 8.7 at 150 kpc and r t /r s = 6.3 at 100 kpc). Due to the 



large ratio, very little matter would be stripped from the satellite 
even after a long time period in a stable orbit. This was confirmed 
by the results quoted in Table[3] where the 150 kpc orbits for the NR 
case showed a 5 per cent change in the satellite's w max - The mass 
removal case without adiabatic expansion (PR) showed a 7 per cent 
change. 

Even in distant orbits, the combined effects of including 
baryons in simulations, namely mass removal with subsequent adi- 
abatic expansion and the inclusion of a baryonic disc, doubled the 
reduction in i> ma x from 5 to 10 per cent (for 150 kpc orbit) and 1 1 
to 22 per cent (for 100 kpc orbit), as given in Table [3] Mass loss 
was even more significant than the change in u max , with losses of 
60 per cent in a 150 kpc orbit and 80 per cent in a 100 kpc orbit for 
the ER case. This shows that the evolutionary history of galaxies in 
pure DM simulations may produce incorrect masses of subhaloes 
at redshift zero, even if they do not interact closely with their hosts. 

Closer circular orbits at 70 and 50 kpc were run to examine 
more extreme tidal stripping. Although the haloes had different 
r t /r s values, all underwent significant mass loss after 5 Gyrs. The 
most dramatic of these were the ER and IR 50 kpc orbits. Fig.|4]dis- 
plays the significant tidal stripping satellites experience when orbit- 
ing close to their host, by examining their circular velocity profiles. 
The ER case is shown every Gyr throughout its 50 kpc circular or- 
bit, and the initial and final IR profiles are shown for comparison. 
This allows one to trace the evolution of the entire circular velocity 
profile for this specific case. 

The v msx and r max of the ER case decreased with time until 
the final satellite has w max = 10 km s _1 , with a significant fraction 
of bound mass existing as tidal arms and tails. At 2 Gyrs, the v max 
has decreased by 40 per cent, converting a 51 km s _1 subhalo to 
a much smaller 30 km s _1 satellite. The initial and final IR veloc- 
ity profiles are also shown in Fig. |4]for comparison. Although the 
ER and IR cases began with similar circular velocity and density 
profiles, the final IR case is slightly more disrupted than the ER 
subhalo, which retained a small satellite core. 
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Figure 6. Mass evolution and fractional change in maximum circular velocity (« m ax/fmax,i) for the pre-halo removal (PR), which is on the left, and 
exponential removal (ER) satellite, which is on the right. The four elliptical orbits: l:2e (150 kpc to 70 kpc), l:3e (150 kpc to 50 kpc), l:4e (120 kpc to 30 kpc), 
and l:5e (150 kpc to 30 kpc) are shown. Mass continually decreases throughout the 6 Gyrs of evolution. Pericentre passages occur at the steep drops in both 
mass and timax- The l:4e and l:5e orbits have the same pericentre distance of 30 kpc, with the difference being their apocentres. The l:5e orbit experiences 
less mass and D max loss at a given time because of the longer time to complete one orbit. If this time-scale is taken into account, the same pericentre distances 
create very similar evolutionary histories. 



Observational data points of the half-light velocities for the 
classical Local Group dSphs are shown in Fig.|4]as well. The initial 
ER profile is too massive to be consistent with the majority of Lo- 
cal Group dSphs but is consistent at 2a with Draco, the data point 
just above it, and Leo II, the data point with the smallest radius 
value. After 1 Gyr of evolution, the ER subhalo is consistent with 
four observational data points, and after 3 Gyrs, two more observa- 
tions match the ER satellite's circular velocity profile. At the end 
of 5 Gyrs of evolution, the satellite is capable of matching all ob- 
served Local Group dSphs at some point in its evolutionary history. 
The gravitational effects due to the expulsion of baryons have con- 
verted a massive satellite, incapable of matching any of the Local 
Group dSphs, into a dwarf that matches a given observed dSph at a 
fixed time in its evolution. 

The ER and IR satellites have a r t /r s = 2.2 for the 50 kpc 
orbit. As seen in the evolution of the circular velocity profiles, 
satellites with this low tidal to scale length value are capable of 
undergoing extreme disruption and possibly becoming completely 
destroyed. It is within the 2 < r t /r s < 4 region that tidal stripping 
may become significant depending on pericentre passage, number 



of orbits completed, and central density of the satellite. This will be 
further examined in our elliptical orbits. 

As a comparison between the satellites run in circular orbits, 
Fig.[5]shows the circular velocity profiles of the NR (left panel) and 
the ER cases (right panel). This figure also includes observational 
data points of the half- light velocities for the classic al Local Group 
dSphs and lines from iBovlan-Kolchin et alj J2012h Fig. 3 for the 
Aquarius E halo. We plot both the isolated NR and ER haloes over 
these data and the results of the 150 kpc, 100 kpc, 70 kpc, and 
50 kpc circular orbits at 5 Gyrs. 

The initial NR subhalo was created in order to match cosmo- 
logical simulations, without baryon removal and without a MW 
disc component. The circular velocity profile for the NR halo is 
consistent with the Aquarius satellites in that it does not match any 
of the observed MW satellites. The initial ER satellite is in better 
agreement because it includes mass removal even before it evolves 
on a circular orbit. This demonstrates how baryonic removal alone 
brings the subhaloes into agreement with MW satellites. Creating 
an initial NR dwarf with a more typical, low total mass value, would 
provide even better agreement between the initial ER dwarf and ob- 
served subhaloes. 
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Figure 7. Shown is the evolution of circular velocity at 500 pc for the expo- 
nential removal (ER) case. Each elliptical orbit (l:2e, l:3e, l:4e, and l:5e) 
is shown with a different linestyle. The velocity at this set radius decreases 
dramatically as the satellite approaches pericentre, similar to how « ma x 
evolved in Fig. [6] 

After undergoing circular orbits for 5 Gyrs, the NR halo re- 
mained inconsistent with observations for the 150 kpc, 100 kpc, 
and 70 kpc orbits. This i s in agreement with the results of 
iBovlan-Kolchin et alj fioH) . who found many subhaloes that were 
too dense and had circular velocities that were too large to match 
the MW satellite population. Only with an extremely close circular 
orbit (50 kpc), does the NR satellites agree with Draco. 

The ER satellite included both mass removal and a MW disc 
component. The results for this halo are in much better agreement 
with the entire population of Local Group dSphs. As seen in Fig.[5] 
the 100 kpc orbit is in agreement with four of the observed dSphs. 
When the satellite is placed in a circular orbit with a 50 to 70 kpc 
radius, it is able to match the four dSphs with the smallest circular 
velocity profiles. The importance of including the MW disc is seen 
from these orbits as significant tidal stripping reduces their mass 
and maximum circular velocity values. 



5.2.2 Elliptical orbits 

Results from our circular orbits motivated us to create elliptical or- 
bits that have pericentres within the interesting range of truncation 
radii from 1 < r\/r B < 5. Our range in dwarf masses, in combi- 
nation with a desire to create realistic ellipticities and pericentres, 
resulted in four orbits with apocentre-pericentres of 150-70 kpc, 
150-50 kpc, 120-30 kpc, and 150-30 kpc. 

As the satellites began their orbits, they immediately started to 
lose mass under the tidal influence of the MW. At their first peri- 
centre passage, most had already lost over half their mass. The sec- 
ond orbit continued to induce strong mass loss, and by their second 
pericentre passage, some satellites had only one fifth of their mass 
remaining. Both the mass and maximum circular velocity evolution 
for the PR and ER satellites over 6 Gyrs are shown in Fig.|6]for all 
different orbits. The maximum circular velocity plot is shown as a 
fraction of the original maximum circular velocity value, while the 
mass is plotted in solar masses. 

Each "step" in the circular velocity plot marks a pericentre 
passage, which also aligns with the "bounce" or change in slope 
of the mass evolution. However, mass is lost continuously through- 



out the orbit, as opposed to v max , which only changes significantly 
at pericentre passage. For the most distant orbit (l:2e), the mass 
shows a smooth decrease, but the closer orbits have increased mass 
loss and change in circular velocity. At pericentre, the mass tends 
to increase slightly before decreasing once again. This is due to 
some of the mass in the tidal features temporarily becoming bound 
to the satellite. We found v mBX to be a more robust parameter for 
describing the transformation a satellite has undergone throughout 
its orbit. 

Note that the l:4e (120-30 kpc) and l:5e (150-30 kpc) lines 
occur in reverse order on the plot, with l:5e undergoing less mass 
and Umax loss than the l:4e case. Although this may seem strange 
at first glance, it is due to the pericentres and apocentres chosen 
for the orbits. The 1 :2e, 1 :3e, and 1 :5e orbits all have apocentres of 
150 kpc. However, the l:4e orbit has an apocentre of 120 kpc and 
the same pericentre as the 1 :5e orbit, resulting in a shorter orbital 
time-scale. This means the 1 :5e orbit has less mass loss at any given 
time. What is most interesting is that after the same number of peri- 
centre approaches, the change in mass and w max is approximately 
the same for the 1 :4e and 1 :5e orbits. This shows that the pericen- 
tre distances, and consequently the tidal radii, are good indicators 
of total tidal stripping. 

A more accessible parameter for observers is a measurement 
of the circular velocity within the central region of the galaxy. To 
this effect, we provide the evolution of circular velocity at 500 pc 
in Fig. [7] Here, we see similar evolution as for w max but with final 
values in the range of 14 to 19 km s -1 . 

When comparing the elliptic al orbit profiles to the Loca l 
Group data points and the work of IBovlan-Kolchin et alj d20 12h . 
as done in Fig. [5] with circular orbits, we find similar trends. The 
NR halo is not compatible with the Local Group dwarfs in ellipti- 
cal orbits. The PR halo with l:2e and l:3e orbits are only consistent 
with Draco, and the 1 :4e and 1 :5e orbits are consistent with two ad- 
ditional Local Group dwarfs within 2a. 

The initial circular velocity profile of the PR halo was lower 
than the NR case but only matched the initial ER halo profile after 
a l:4e orbit around the MW. We find that simply including a disc 
component in the MW is not sufficient to create agreement between 
our dwarf and the Local Group dSph population. Both the PR and 
ER haloes have the same total mass and evolve on the same orbits, 
providing evidence that the internal structure of the dwarfs drives 
the difference between their final profiles. 

The ER halo is able to match all but the lowest two Local 
Group dwarfs with the elliptical orbits inducing significant tidal 
stripping. The results of the ER l:2e and l:3e orbits are similar 
in strength, lying between the 100 kpc and 70 kpc circular orbits, 
for the same halo. The l:4e and l:5e orbits are stronger than the 
70 kpc circular orbit but not as dramatic as the 50 kpc orbit, with 
fcirc(Ll kpc) = 20 kms -1 for the l:4e orbit. This agreement 
with the Local Group dwarfs once again shows the importance of 
including our two baryonic effects within simulations in order to 
match observations. 

Next, we investigated how changing the ellipticity and peri- 
centre distances affected mass removal, based on r t /r s . In Fig. [8] 
we show the fractional change in mass and v max after 5 Gyrs in 
relation to the r t /r s of each satellite. We calculated the truncation 
radius of our subhaloes at pericentre according to equations dl lb 
and J 1 2b . Circular orbits have radii of 150 kpc, 100 kpc, 70 kpc, 
and 50 kpc, and elliptical orbits have pericentres of 70 kpc, 50 kpc, 
and 30 kpc. 

To evaluate the truncation radii values, we input the analytic 
form (disc plus NFW) of the MW mass, M(R), and pericentre dis- 
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tance, R, into equation i ll lb . In order to obtain a correct value for r t , 
no analytic form of the satellites' mass profile, m(r), was used. In- 
stead, numerical results from the isolated equilibrium haloes were 
input. The tidal radius, r t id c , was then solved for each satellite and 
pericentre combination. We also solved for the resonant radius us- 
ing this method. The minimum of r Tes and r t id c was taken as the 
truncation radius, r t . This was divided by the satellites' r s to create 
the value r t /r s . 

Now that we have defined r t /V s for each satellite case with 
a different pericentre, we examined how the tidal radius impacts 
the mass loss after 5 Gyrs of evolution. In Fig. [8] we show satellite 
cases in different colours and symbols for different orbits, as listed 
within the figure. The instantaneous removal case is not included 
on this plot because the mass and u max loss versus r t /r s values are 
almost identical to those seen in the exponential removal case. 

It is apparent that for satellites with pericentres that are large 
enough to reach r t /V s > 4 the tidal stripping is minimal. As ex- 
pected, over half of a satellite's mass is lost, but the v max does 
not suffer a similar reduction. The fraction of i> ma x remains above 
0.8 for these cases. For example, the shift in r t /r s between the 
NR 150c and ER 150c points is due to the adiabatic expansion of 
the mass removal haloes, which increases their r s values. This gives 
smaller r t /V s values, while only slightly lowering their i) max frac- 
tion. 

However, for satellites on orbits that reach r t /V s ~ 4 and 
below, there is a broad range of possible mass loss scenarios. In 
this regime, a satellite's mass can range from 0.3 to 0.05 times the 
original total mass, and the circular velocity can drop to 0.8 to 0.2 
its original value. This is a dramatic decrease, indicating that tidal 
radius over scale radius at pericentre is an important parameter for 
predicting mass loss. 

The rt /r s values of the PR haloes are significantly larger than 
those of the ER haloes, as denoted by different colours with the 
same symbols. The internal structure of the halo makes the PR 
cases resistant against tidal stripping so that the satellite does not 
experience as large a reduction in w max or mass. 

The separation between the circular and elliptical orbits is due 
to the length of time that the satellites sit near the primary. As cir- 
cular orbits are always at the same distance, they are able to reduce 
their mass and u max more, over 5 Gyrs of evolution, than an el- 
liptical orbit with the same pericentre. Points with the same r t /r s 
value indicate that they have the same pericentre. This makes for an 
easy comparison between points that lie on vertical lines. For ex- 
ample, the 70 kpc circular orbit (70c) and the 150-70 kpc elliptical 
orbit (l:2e) have the same r t /r s for each satellite case of a single 
colour. The vertical separation between these pairs of points is due 
the length of time they are close to the primary. 

Although satellites on circular orbits undergo a larger number 
orbits around the MW in the same amount of time as the elliptical 
orbits, the tidal forcing they feel is constant throughout all time. 
This means that tidal truncation is the main method of mass re- 
moval, without tidal shocks. We have included our tidal stripping 
model, as described in Section |3~4l that mimics tidal truncation. 
Unbound mass is iteratively removed from the satellite, and the 
tidal radius is recalculated with each iteration. The results of this 
model with 10 and 30 iterations are plotted in Fig. [8] Our simula- 
tion results typically fall between these two lines, indicating that 
the number of pericentre passages, as modelled by the number of 
iterations, can account for the spread between points with the same 
r t /r s value. 
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Figure 8. Results of tidal stripping for the mass and Umax of a satellite af- 
ter 5 Gyrs. Symbols indicate different orbits with open symbols for circular 
orbits and filled symbols for elliptical orbits. The colours indicate different 
satellite cases: red is for exponential removal, blue is for pre-halo formation 
removal, and black is for no removal. Different satellite cases (colours) with 
the same orbit (symbol) always occur from left to right as ER, PR, and then 
NR, with increasing values of rt / r B . The fraction of mass retained increases 
linearly with increasing r t / r s so that the NR case always has the largest fi- 
nal mass fraction. When rt / r s is below 4, tidal stripping is very important, 
as seen by the significant mass loss and decrease in f max f / v max i ■ The 
solid black lines display our model for tidal stripping by iteratively remov- 
ing all unbound particles and recalculating the tidal radius. The results of 
10 and 30 iterations are the top and bottom solid lines, respectively. 



6 CONCLUSIONS 

One of the goals of this paper is to resolve some of the small-scale 
problems of the standard ACDM model, namely the apparent over- 
abundance of satellites with 30 k m s -1 < v c i TC < 50 kmr 1 
jBovlan-Kolchin et alj20lll 20121) . In line with other recent results 
of cosmological simulations of galaxy formation dZolotov et alj 
120121 : iBrooks & Zolotovll2012l ; IBrooks et alj|2012l) , we argue that 
the solution of the problem is in the gas physics of forming dwarf 
galaxies. By posing the question what are the effects of baryons on 
dwarf galaxies, we hope to examine just one aspect of the increas- 
ingly complicated fields of galaxy formation and evolution. 

Most of the existing theoretical predictions and constraints on 
the abundance and structure of small satellites ("substructures") are 
made using iV-body simulations. Because observed dwarf galaxies 
in the MW are heavily dominated by the dark matter, one expects 
that there is little impact of baryonic physics on the structure of 
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dwarf satellites. However, this appears not to be true; the correc- 
tions due to baryons are significantly larger than one would naively 
expect. 

Baryon removal may reduce the maximum circular veloci- 
ties of satellites enough to explain the "too big to fail" problem. 
The central mass distribution of a dwarf may be substantially al- 
tered both by gas removal at infall and by tidal stripping due to the 
host galaxy. iV-body simulations underestimate stripping because 
(i) baryons are not removed from the dwarf and (ii) the mass of the 
host is not enhanced by infalling baryons. Accounting for these ef- 
fects can alter the circular velocity profiles of the subhaloes, such 
that they are "small enough to succeed" in explaining the observed 
structure of the Milky Way dSphs. 

We used a number of methods to investigate how baryonic 
physics impacts the structure of a dwarf spheroidal. We ran hy- 
drodynamical cosmological simulations, high resolution TV-body 
simulations, and an adiabatic expansion model to examine the ef- 
fects of baryon removal from dwarf galaxies. The impact of tidal 
stripping on dwarf galaxies was explored using a series of high res- 
olution TV-body simulations and two analytical models. 

Our two high resolution cosmological simulations were run 
with and without gas removal. The gas removal case was created to 
reproduce the heating and removal of gas as a dlrr falls into a host 
galaxy. This dwarf expanded adiabatically in response to the mass 
loss. Baryon removal from the TV-body simulations also resulted 
in adiabatic expansion for our exponential and instantaneous mass 
removal cases. The rate of removal of the gas does not affect the 
final results. In other words, both fast and slow removal provide the 
same amount of expansion. 

The gas re mova l results were modelled with the 
iBlumenfhal et all Jl986h prescription for adiabatic expansion. 
Mass removal causes dwarfs to both expand and lower their 
circular velocities. The magnitude of the expansion depends on 
the fraction of mass removed and on the slope of the density. In 
the case when the same fraction of mass is removed independent 
of radius, the shape of the density profile does not change. The 
maximum decline in the density, according to adiabatic expansion 
model, is equal to (1 — /t,) 4 , where fy is the cosmological fraction 
of baryons. This decline is expected to happen for galaxies with 
constant density cores. For galaxies with cusps, the decline is 
smaller. For the NFW = — 1 cusp, the decline is still substantial 
with pi/ Pi = (1 - /b) 3 . If we assume that the mass removed is 
the cosmological average fraction, then pf/pi ~ 0.5. 

We tested a series of circular and elliptical orbits for each 
dwarf with high resolution TV-body simulations. Without baryon 
removal, tidal stripping duri ng an orbit around a pure-NF W host 
creates profiles similar to the lBovlan-Kolchin et al.l f2012l) results. 
Our no mass removal dwarf does not agree with the half-light ve- 
locities of the Local Group dSph population after 5 Gyrs of evolu- 
tion in circular or elliptical orbits. 

The combination of tidal evolution around a MW with both 
a NFW and disc component and baryonic removal with adiabatic 
expansion was tested with the exponential and instantaneous mass 
removal dwarfs. Both experienced large amounts of mass removal 
and a significant reduction in v max . This brought them into agree- 
ment with every Local Group dSph over the course of at least one 
of the elliptical or circular orbits. 

Tidal evolution alone was also examined in our pre-halo for- 
mation mass removal dwarf. This subhalo had the same mass as 
our other baryon removal cases but without adiabatic expansion. 
It experienced weaker tidal stripping due to its steep inner density 
profile, making it compatible only with the five most massive Lo- 



cal Group dSph population after evolution in elliptical and circular 
orbits. 

In order to quantify how the inner regions of the dwarf are im- 
pacted due to tidal stripping, we created two analytical models. The 
first calculated how a dwarf's tidal radius changes depending on the 
mass of the dwarf, the mass of the primary, the distance from the 
main halo, and the inclusion or absence of baryons in the central re- 
gion of the host. We found that r% /r s , the tidal radius over the scale 
radius of the satellite, is able to predict the strength of tidal strip- 
ping. When rt/r a is greater than 4, there is no significant evolution 
in t) mal . However, when r t / r s is less than two, the central region of 
the dwarf is subject to tidal forces and may be significantly altered. 

The second analytical model for tidal stripping iteratively re- 
moved unbound particles, estimated the correct adiabatic expansion 
the halo would undergo from this amount of removal, and then re- 
calculated the tidal radius. This model was able to reproduce our 
TV-body simulations in the context of reduction in D max and mass 
as a function of r t /r s and time in orbit. We find that each iteration 
continues to remove some fraction of the particles without conver- 
gence on a value. 

Our main conclusions are as follows: 

(i) Our simulations show that without b aryon removal or tidal 
strippi ng, dwarfs match the results from the iBovlan-Kolchin et alj 

analysis of the Aquarius simulations. 

(ii) Baryon removal with adiabatic expansion is very effective, 
bringing initially extremely massive dwarfs into agreement with 
observations before any orbital evolution. 

(iii) Including the effects of tidal stripping without expansion 
due to baryon removal is insufficient to reproduce the entire ob- 
served dSph population for the case of an initially massive dwarf. 

(iv) The inclusion of both baryon removal with adiabatic expan- 
sion and tidal stripping is extremely powerful, capable of trans- 
forming an initially incompatible dwarf into one that agrees with 
the entire Local Group dSph population. 

(v) Dwarfs undergo continual mass removal and reduction in 
Umax at each pericentre passage. These values do not converge even 
after 6 Gyrs of evolution. 

(vi) Tidal radius over scale radius, r t /r s , is an important param- 
eter for understanding the potential mass and « max loss of a dwarf. 
If r t /r s is less than four, satellites may become substantially dis- 
rupted, depending on the number of times they orbit their host. 
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